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Theory  of  Propagation  of  Electromagnetic 
Wgy®8  in  Space-Time  Varying  Media 


5.  INTRODUCTION 

A  difiicult;  problem  to  consider,  conceptually,  is  wave  propagation  in  media 
which  vary  with  both  position  and  time.  In  such  media  our  standard  concepts  of 
frequency,  wavenumber,  and  group  velocity  no  longer  apply.  That  is,  we  find  that 
frequency  can  only  be  defined  as  the  time  derivative  of  the  phase  function,  and  that 
the  quantity  so  defined  is  a  function  of  both  the  position  and  the  time  at  which  the 
observation  is  made.  We  also  find  that  the  group  velocity  no  longer  retains  its 
conventional  meaning.  In  spatiaily-homogeneous,  time -invariant  media  the  group 
velocity  is  interp  *eted  as  the  velocity  at  which  wave  packets  centered  around  some 
wavevecto?  kj,  propagate  (Jeffreys  and  Jeffreys,  1962),  In  space-time  varying 
media  this  is  no  longer  true.  La  fact,  as  we  shall  see,  values  of  u  and  k  do  not 
propagate  with  the  group  velocity  Vj,w  ;  rather  different  quantities  (which  a~e 
functions  of  w  and  k)  are  propagated  at  this  velocity.  The  same  conclusions  hold 
true  for  energy  flow  (that  Is,  the  energy  flux  does  not  propagate  with  the  group 
velocity). 

In  this  paper  we  will  study  the  propagation  of  electromagnetic  waves  in  loss¬ 
less  media  which  vary  slowly  with  position  and  time.  We  will  therefore  employ 
the  four- dimensional  WKB  (’'Ver.tr.el-Kr  smrrs-Brillouin)  method.  The  WKB  method 
was  first  applied  ir.  three  dimensions  by  Sommerfeld  and  Runge  (1931).  That  is. 
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in  media  in  which  the  properties  depended  on  position  x  but  not  on  time, 

Sommerfeld  and  Runge  considered  solutions  of  the  form  exp  [iwt  -  Jk  •  dxj. 

Because  the  phase  function  must  be  uniquely  defined,  Sommerfeld  and  Runge  then 
concluded  that  7  X  k  =  0  was  required.  This  is  the  original  version  of  the 
Sommerfeld -Runge  Law.  This  result  can  be  extended  to  the  four-dimensional  case, 
as  has  been  done  by  Whitham  (1960)  and  Poeverlein  (1962).  That  is,  in  space-time 
varying  media  the  existence  of  a  uniquely  defined  wave  function,  exp  [i f(u dt  -  k •  dx)J 
requires  that  7w  +  d  k/dt  =  0.  This  is  the  four- dimensional  Sommerfeld -Runge  Law. 

In  the  present  paper  we  will  first  give  an  elementary  denv  Uion  of  the  four¬ 
dimensional  Sommerfeld-Runge  Law.  We  will  then  examine  its  implications,  and 
finally  indicate  its  use  in  studying  electromagnetic  wave  propagation  in  isotropic, 
lossless  media  which  vary  slowly  with  position  and  time. 


2.  GENERAL  THEORY 


2. 1  Derivation  of  the  Generalized  Sommerfeld-Runge  Law 

To  introduce  the  Sommerfeld-Runge  Law  we  begin  with  the  pair  of  equations 
for  the  electric  field  E 


VxVxe  =  -u 

—  ro 


d2(«  E) 
at  2  ’ 


(la) 


7  ’  («E)  =  0, 


(lb) 


where  the  permittivity  « (x,  t)  is  a  function  of  both  position  and  time  but  the  mag¬ 
netic  permeability  is  equal  to  that  of  vacuum.  Equations  (la)  and  (lb)  are  appropri- 
ate  fer  propagation  in  a  lossless  space -time  varying  die’ectric.  We  will  later 
generalize  our  results  to  include  dispersive  media,  such  as  a  plasma.  Let  us 
assume  a  solution  for  E  of  the  form 


E(x,  t)  =  eo(x.  t)  e  ^ l)  ,  (2) 

where  eQ  varies  slowly  with  position  and  time  compared  with  <£{x,  t).  Using  (2)  in 
(lb)  we  find  that,  if  i  j  7«  |  <<  |  7t/>|  we  get 

7<#>-eo  =  0.  (3) 

If  we  use  (2)  and  (3)  in  (la)  we  then  find  that  if  the  dielectric  varies  slowly  with 
position  and  time,  the  function  4>  satisfies 


=  /xQ  «(x,  t) 


where  the  conditions  that  the  medium  be  considered  slowly  varying  are 
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(4) 


|  V*  I2  »  V2<f>  , 

(5a; 

[dtj  dt2* 

(5b) 

V 

V 

*1“ 

(5c) 

l^i  »|  |v<|. 

(5d) 

The  solution  '  Eq.  (4)  can  be  written  quite  generally  as 

<p  ■  i  J{  k  •  dx  -  w  dt)  , 

L  “ 

(6) 

where  L  is  a  line  integral  in  four  space  between  some  initial  point  (x^  tQ)  and  the 
observation  point  (x,  t).  We  can  identify  :o  and  k  in  (6)  by  realizing  that  since  the 
phase  function  <f>  must  be  a  single  valued  function,  the  line  inir-gral  in  (6)  must  be 
independent  of  the  path  joining  (x^,  tQ)  to  (x,  t).  Therefore,  the  integrand  in  (6) 
must  be  a  perfect  differential  and  we  have 

ik  -  ( V<£  )t  ,  (7a) 


Using  these  expressions  in  Eq.  (4)  gives  the  dispersion  relation  between  w  and  k  as 

c2k2=«r«2,  (8) 

2 

where  k  =  ki  k  and  «r  is  the  relative  permittivity  of  the  dielectric.  Equations  (7a) 
and  (7b)  imply  a  relation  between  w  and  k,  in  addition  to  that  given  by  Eq.  (8). 
Differentiating  (7a)  with  .  aspect  to  t,  taking  the  gradient  of  (7b),  and  then  adding 
the  resulting  equations  gives: 


Equation  (8/  is  the  four-dimensional  version  of  the  Sommerfeld-Runge  Law,  This 
result  can  alternately  be  derived  by  employing  arguments  on  wave  conservation,  as 
has  been  done  by  Lighthill  and  Whitham  (1955).  Equation  (9)  has  also  been  dis  - 
<  ^ssed  by  Whitham  (1960),  Poeverlein  (1962),  and  Landau  and  Lifschitz  (1958), 

2.2  Discussion  of  the  Properties  of  Eq.  (9) 

2.  2. 1  GENERALIZED  GROUP  VELOCITY 

To  discuss  the  properties  of  Eq.  (9)  in  a  medium  which  varies  with  beta  posi¬ 
tion  and  time,  let  us  first  write  Eq.  (8)  in  the  general  form 

w  =  W  (k,  x,  t)  (10a) 


or 


k  =  K (a),  x,  t)  .  (10b) 

Equations  (10a)  and  (10b)  are  valid  for  an  arbitrary  linear  medium,  and  not  just  a 
dielectric.  For  example  in  a  lossless  plasma 

k  =  c"1  •  -v/w2*  Wp  (x,  t) 

where  «»0  is  the  electron  plasma  frequency.  If  (10a)  is  substituted  into  (9)  we 
obtah. 

!f+(V.V)k  =  -(VW)M  (11) 

where 

i"<7kW)‘  (It).  «»> 


The  quantity  V  in.  (12)  .tar,  be  interpreted  as  a  generalized  group  velocity  as  we 
shall  see  in  the  following  discussion.  UKon  defining  d/dt  =  d/dt  +  V*V  Eq.  (11) 
can  be  rewritten  as 


5 


dk 

— gx 

dt 


-  (VW) 


k,  t  * 


(13) 


From  (13)  we  see  that  if  W  does  not  depend  explicitly  on  position,  then  dk/dt  =  0. 
This  means  that  if  one  moves  along  the  ray  with  the  velocity  V,  ho  will  observe 
constant  values  of  k  (that  is,  the  observer  moving  with  V,  measures  constant 
values  of  wavelength).  Therefore,  if  W  is  a  function  of  k  and  t,  but  does  not 
depend  explicitly  on  x,  then  k  is  a  constant  of  the  motion,  so  that  wave  packets* 
sharply  centered  around  some  wavenumber  kg  will  be  propagated  with  V  =  dW/dk 
evaluated  at  k  =  k  . 

—  — -O 

To  examine  the  other  limit  when  the  properties  of  the  medium  vary  with  posi¬ 
tion,  but  not  with  time,  it  is  convenient  to  use  (10b)  in  (9).  For  this  case  we  obtain 
(assuming  the  medium  is  isotropic) 


t V • Vw 
at  — 


x 


(14) 


From  Eq.  (14)  it  is  clear  that  if  K  depends  on  w  and  x,  but  does  not  explicitly 
depend  on  time,  then 


dw 

dt 


0. 


(15) 


This  means  that  the  observer  moving  along  the  ray  with  the  velocity  V  will  observe 
constant  values  of  w  (that  is,  constant  wave  period;,  Snell's  Law  follows  immedi¬ 
ately  from  Eq.  (15).  Therefore  wave  packets  sharply  centered  about  some  fre¬ 
quency  uQ  will  propagate  with  the  group  velocity  V  (evaluated  at  u>  =  wQ). 

For  the  general  case,  when  the  properties  of  the  medium  depend  on  both  posi¬ 
tion  and  time  neither  o>  nor  k  will  be  invariant  as  one  moves  along  the  ray  with  the 
group  velocity  V.  (In  Section  2. 3,  however,  we  do  show  that  there  are  velocities 
V(k)  \  V  and  yfa  ^  H  V  with  which  values  of  k  and  hi  are  propagated  in  space-time 
varying  media. )  In  this  case  other  quantities  will  be  invariants  of  the  motion. 
Consider  each  scalar  component  of  Eq.  (11).  We  have 


fUj. 

dt 


+  (V  •  V 


)kl  = 


(16) 


*In  a  spatially  homogeneous,  time -varying  medium,  pulses  can  be  represented 
as  Jdk  A'k)  exp[,(k*  x  -  Jw  (k„  t)  dt)  J  provided  A(k)  has  significant  amplitude 
only  at  values  of  k  for  which  Eq.  (5)  is  satisfied. 
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To  solve  Eq.  (16)  we  consider  the  subsidiary  set 


Let  us  denote  the  particular  integrals  of  (17)  by  f(k^,x, t)  =  Cjj  g(k^,x,  t)  =  C,,; 
Mk^.x.t)  =  C3;  i^(k%,x,  t)  =  .  Then  it  can  be  shown  that  (Sneddon,  1957) 


( s )  The  general  solution  of  (16)  is  given  by 

Cj  •  <t>(C2,  C3,  C4)  .  (18) 

where  i>j  an  arbitrary  function,  determined  by  the  boundary  conditions  imposed. 

(b)  The  invariants  of  motion,  for  an  observer  moving  with  the  velocity  V  given 
in  Eq.  (12),  are  Cj,  C 2«  Cg,  C^.  That  is,  each  Cg  satisfies 

dC  i  ,  i 

■a r  Mir+x*7)  cs  =  °-  (19> 

Therefore,  for  media  which  vary  with  both  position  and  time  the  generalized 
group  velocity  defined  in  Eq.  (12)  is  the  velocity  with  which  the  quantities  Cs  are 
propagated.  It  is  only  in  the  limit  of  spatially  homogeneous,  time-invariant  media 
that  some  of  the  constants  C„  can  be  identified  with  w  and  k . 

In  Section  8. 1  we  will  calculate  the  constants  of  motion  for  several  examples 
of  media  which  vary  with  position  and  time . 

2.  2.  2  THE  ANGLE  BETWEEN  V  AND  k 

In  this  section  we  will  demonstrate  that  in  isotropic  lossless  media  the  normal, 
£  =  (k/k),  to  the  phase  surface  4>  (x,  t)  lies  in  the  same  direction  as  the  group 
velocity  V.  We  will  alsr  calculate  the  angle  between  fr  and  V_  for  the  case  of  a 
simple  isotropy. 

Let  us  consider  UOb)  for  the  case  when  the  dispersion  relation  can  be  written 
as 


k  =  K(  a.,  g.x,  t),  (20) 

where  d  is  the  angle  that  the  vector  k  makes  with  the  z  axis.  (This  corresponds 
to  the  case  when  the  dispersion  relation  is  given  by  the  Appleton- Hartree  formula. ) 


Using  Eq.  (20)  in  Eq.  (12),  we  have  upon  differentiating  implicitly  that  the  group 
velocity  is 


i-HH)M(flf)]- 


where  k  is  a  unit  vector  along  k,  and  £  is  the  unit  8  vector  in  a  spherical  coordin¬ 
ate  system.  Taking  the  cross  product  of  k  with  (21)  gives 


**.v  I  a(fr)  (w)- 


om  (22)  we  see  that  if  the  medium  is  isotropic  and  lossless  so  that  K  is  inde¬ 
pendent  of  8  in  Eq.  (20),  then  k  x  V  =  0  so  that  k  and  V  are  in  the  same  direction. 
The  dot  product  of  k  with  (21)  gives  k  •  V  =  k  (  dco  /  d  k) .  From  this  result  and  (22) 
we  can  compute  the  angle  y  between  k  and  V .  We  have 


*  I  1**21 
tan  y  |  -  k .  y 


1  dK  1  dk_ 
k  39  '  k  69  * 


Equation  (23)  is  a  well-known  result  from  the  study  of  whistlers  (Kelt  and  Haskell, 
1965;  Kelso,  1964),  We  have  shown  here  that  the  result  also  applies  to  media  which 
vary  slowly  with  both  position  and  time,  provided  we  understand  that  V  is  not  the 
velocity  with  which  values  of  k  and  u>  are  propagated.  One  final  result  of  interest 
is  to  use  (23)  and  (21)  to  calculate  the  magnitude  of  the  group  velocity.  We  get 


m  «/ 


2.  2.  3  THE  EQUATION  OF  MOTION  OF  £ 

In  this  section  we  shall  demonstrate  that  the  classical  equation  presented  by 
Landau  and  Lifschitz  (1959)  for  the  motion  of  the  normal  to  <j>  (x,  t)  can  be  general¬ 
ized  to  include  space- time  varying  media.  For  isotropic  media  we  write,  using 
k  =  k  ic  in  Eq.  (1 1), 


We  next  use  the  dispersion  relation  of  Eq.  (10b)  in  the  second  term  on  the  left  hand 
side  of  (25).  We  obtain,  after  regrouping  terms 
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where  V  =  |  V  |  =  {dK/du  J*1 .  Now  substitute  Eq.  (14)  for  dw  /dt  into  Eq.  (26). 

The  result  is 

k^*V«.  (7K)„_,.-(VW)k_t.  (27) 

Equation  (27)  is  valid  for  an  arbitrary  lossless  isotropic  medium.  To  obtain  the 
analog  of  the  ray-normal  equation  of  Landau  and  Lifschitz  (1959',  we  next  special¬ 
ize  (27)  to  dielectrics.  For  this  case  K  =  w/v  where  v  =  V  =  phase  velocity.  In 
this  limit  Eq.  (27)  becomes 

^  -  £<£  •  Vv)  =  -  Vv(x,  t)  .  (28) 

Equation  (28)  determines  the  motion  of  the  ray  normal  in  dielectrics  which  vary 
slowly  with  both  position  and  time.  As  demonstrated  by  Landau  and  Lifschitz,  who 
obtained  the  same  equation  for  the  case  when  v  varies  with  position  only,  Eq.  (28) 
predicts  a  bending  of  the  rays  toward  the  region  where  v  is  smaller. 

We  also  note  from  Eqs.  (27)  and  (28)  that  if  v  is  independent  of  position,  but 
does  depend  on  time,  we  have 

=  0.  (29) 


Therefore,  as  expected,  the  ray  does  not  change  its  direction  of  propagation  in 
media  that  vary  only  with  time. 

2.  2.  4  TEMPORAL  DISCONTINUITIES 

It  is  often  desirable  to  know  the  behavior  of  w  and  k  when  the  properties  of  the 
medium  are  suddenly  altered.  For  example,  suppose  we  have  a  dielectric  in  which 
«  =  <fj  (x)  for  t  <  tj,  and  «  =  SgM  (or  t  >  tj  .  To  study  the  behavior  of  w  an.',  k 
when  temporal  discontinuities  occur  let  us  integrate  Eq.  (9)  from  tj  -  8  to  tj  +  8  . 
(Note  that  Eq.  (9)  is  not  strictly  valid  for  8  =  0  , )  We  obtain 


k  (x,  tj  *-  S  )  -  k  (x,  tj 


V8 


V- 


j  (  Vw  )  dt  . 


(30) 


In  the  limit  at'  8—0,  the  right  hand  side  of  (30)  vanishes  (unless  Vw  has  a  delta 
function  behavior).  Therefore 
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k  (x,  tj  +  8  )  =  k  (x,  tj  -  8  )  . 


(31) 


Since  Ec\  (31)  implies  that  both  the  magnitude  and  direction  of  k  cannot  change 
instantaneously,  then  from  (10b)  we  may  write 

Kj  [w (x,  tj  +8  >,  x,  tj  +  8  ]  =  K2  (  w(x,  tj  -  8  ).  x,  tj  -  8  ]  .  (32) 

For  a  dielectric,  in  which  K  =  ±  u  (x,  t)  ^  «  (x,  t)  , 

Eq„  (32)  yields 


ai  (x,  t  +  8  ) 


tr«i<*> 
"  L  *2 


1/2 


w(x,  t  -  8  )  . 


The  positive  sign  in  Eq.  (33)  is  appropriate  for  the  wave  travelling  along  k,  where¬ 
as  the  negative  sign  is  appropriate  for  the  reflected  component  which  travels  along 
-k.  This  latter  component  is  negligible  in  the  WKB  limit.  In  the  limit  of  a  spatially 
homogeneous  dielectric,  Eq.  (33)  reduces  to  the  previous  result  of  Morgenthaler 
(1958). 


2. 3  Velocity  of  Propagation  of  Value*  of  k  and  u 

One  of  the  points  we  have  made  is  that  in  a  space-time  varying  medium  neither 
w  nor  k  is  an  invariant  as  one  moves  with  the  group  velocity  V  =  Vkw  ,  it  is  possi¬ 
ble,  however,  to  define  a  new  velocity  V<k>  such  that  k  will  be  constant  when  the 
observer  moves  with  this  velocity.  To  obtain  V(k>,  let  us  suppose  that  the  solution 
of  Eq.  (11)  is  given  by  k(x,  t).  Then  we  can  define  a  velocity  V1  through  the 
equation 

(V1  •  V  )  k  =  <VW)^  t .  {34) 

If  this  were  done,  Eq.  (11)  could  be  rewritten  as 
d  k 

~dT+  [(y.+z )  ■  7i!:o.  (o 
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from  which  we  immediately  identify 

y<k*  =  V  +  v'  .  (36) 

Equation  (34)  can  be  solved  for  V  by  standard  matrix  methods.  It  is  interesting  to 
consider  the  one  dimensional  limit  (that  is,  d  /dx d/ dy  =  0).  Then  Eq.  (34) 
becomes 

T' <*•»({?),■  0,1 

which  is  readily  solved  for  V.  Using  this  result,  along  with  Eq.  (12),  in  Eq.  (36) 
then  gives 


(k) 

Therefore,  we  have  demonstrated  it  is  possible  to  define  a  velocity  V  such  that 
the  observer  moving  with  this  velocity  sees  constant  values  of  k.  A  similar  argu¬ 
ment  holds  for  w.  That  is,  we  can  find  a  velocity  v(“)  with  which  values  of  w  are 
propagated  in  space-time  varying  media.  In  general,  ^  will  not  equal  V^', 
except  in  the  limit  of  spatially  homogeneous,  time-invariant  media. 


2.4  Co&ments  on  Energy  Flow 

Since  in  a  space-time  varying  medium  w  and  k  are  not  propagated  with 
V  =  Vkw ,  we  should  not  be  surprised  to  find  that  the  energy  flux  does  not  flow  with 
this  velocity  either.  To  consider  this  problem  let  us  limit  cur  discussion  to  iso¬ 
tropic,  lossless  dielectrics.  We  then  have,  upon  using  Eq.  (2)  in  (1)  and  employ¬ 
ing  the  assumptions  expressed  in  Eq.  (5),  that  the  equation  satisfied  by  is 


k  ( 7  ^  •  e_)  +  e  (  V  •  k)  +  2  (k  •  7  )  e 

— •  — O  — O  — «  — •  — 


=  -  fj. 


+  2w 


+  2to 


(39) 


where  f  =  In  «  .  (Note  that  Eq.  (34)  reduces  to  the  results  of  Section  3. 1.  3  in  Born 

and  Wolf  (1959)  in  the  limit  when  the  medium  does  not  depend  on  time. )  Now  recal- 

$ 

ling  that  k  •  e^  =  0,  taking  the  scalar  product  of  ^  with  Eq.  (39),  and  then  adding 
its  complex  conjugate  to  the  resulting  equation  gives: 
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(40) 

where  X=  £•  x.  If  we  now  substitute  k  =  u  J  fiQt  into  Eq.  (40),  and  use  Eq.  (14) 
to  write  dw/dt  =  -  ^  d'J'/dt,  we  obtain 


JL 

dl 


(41) 


Let  us  define  an  energy  flux  I 


=  1 


2 


Then  in  terms  of  I  Eq.  (41)  becomes 


dl 

Tt 


+  V 


dl  rd± 
H~  ~  I<Jt  • 


(42) 


We  now  consider  the  flux  equation  further  in  the  limit  when  V«  and  k  are  in  the 
same  direction  (say  along  the  z  axis).  From  (42)  we  then  have  that: 

(a)  If  « (z,  t)  is  independent  of  time  then 


dl 

dt 


s  (JL  +  y 

vdt+  v 


0  , 


(43) 


so  that  the  energy  flux  I  remains  constant  as  one  moves  with  the  group  velocity 
-1/2 

V  =  ( « fi  )  .  That  is,  values  of  I  are  propagated  with  the  group  velocity. 

(b)  If  t  (z,  t)  is  independent  of  position,  Eq.  (42)  may  be  written  as 

-5F<‘i>  =  <!t  +  vA)<«i)  =  o*  <44> 


so  that  it  is  not  the  flux,  but  rather  values  of  (*  I)  which  are  propagated  with  the 
group  velocity  when  «  depends  on  time  only. 

(c)  When  the  medium  varies  with  both  position  and  time  neither  I  nor  («I)  are 

1/2 

propagated  with  V.  As  an  example,  suppose  ( nQ*  )  -  /?Q(z)  +  t  /?j(z).  Then  it 

is  readily  shown  that  the  invariant  of  the  motion  is 


I  exp 


/3j <z* )  dz'  . 


(45) 


Therefore,  we  have  demonstrated  that  it  is  only  lor  media  which  are  time  invariant 
that  the  energy  flux  (Poynting  vector)  is  propagated  with  the  group  velocity. 


t*t*r***rctrrfix* 
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For  dielectrics  which  vary  with  position  and  time  it  is  not  the  energy  flux  I,  but 
rather  some  new  quantity  f(I,  x,  t)  which  is  propagated  with  the  generalized  group 
velocity,  given  by  Eq.  (12). 

2.5  Approximate  Solution  of  Gq.  (14)  (Small  Doppler  Approximation) 

In  many  instances  it  is  difficult  to  obtain  exact  solutions  of  Eqs.  (11)  and  (14) 
for  space-time  varying  media.  This  is  especially  true  when  dispersion  is  present. 
However,  when  the  frequency  shift  in  "propagating  through  the  medium  is  small 
compared  with  the  transmitter  frequency,  it  is  possible  to  solve  Eq.  (14)  by  itera¬ 
tion.  Let  us  consider  the  limit  when 


(46) 


where,  as  before  k  =  K(w,  x,  t)  and  ds  =  k  •  dx.  The  integral  is  along  the  ray  path 
in  the  medium  between  the  transmitter  and  the  observer.  For  the  special  case  of  a 
dielectric,  this  requires  that,  in  Edition  to  the  condition  that  the  medium  vary 
slowly,  the  path  length  in  the  medium  cannot  be  so  large  that  J  Kds  »  1.  When 
Eq.  (46)  is  satisfied,  we  may  neglect  the  right  hand  side  of  (14),  and  conclude  that 
oj  is  approximately  a  constant  of  the  motion.  Furthermore,  if  (as  is  usually  the 
case)  we  specify  that  cj  =  wq  for  all  time  at  the  position  of  the  transmitter,  then 
u  a  ojq  for  all  x  and  t  (for  which  Eq.  (46)  is  still  satisfied).  Using  u  =*  u»o  then 
gives  in  Eq.  (4) 

|W>|  =  i-2  n(  «o.  x.  t)  ,  .  (47) 

where  n  (  u,  x,  t)  is  the  index  of  refraction,  defined  by  j  k  |  =  -^  n.  Equation  (47)  is, 
of  course,  the  standard  Eikonal  used  in  ray  optics.  Similarly,  in  Eq.  (6)  we  have 
(substituting  |  k  |  =  3 1  n  (  x,  t) 

riw,  1 

<f>  -  i  [  fn  (  ojq,  x,  t)  tc  •  dx  -  wQt  j  ,  (48) 

where,  as  before,  £  is  the  unit  vector  normal  to  the  phase  surface,  and  is  deter¬ 
mined  by  solving  (47).  Finally,  using  (48)  in  Eq.  (7b)  gives  for  the  first  iteration 
to  the  instantaneous  frequency 

«o  fdn{  V  -X-'  0  , 
w=i  wo  '—J - 5t -  ds 


(49) 


’  *■  ur » P^tB^>^«rirw«^T»»ij*r>  •>  rorow 
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where  ds  =  ft  •  dx  =  path  length  along  the  ray.  Equation  (49)  is  the  result  used  by 
ionospheric  reoearchers  in  studying  the  Doppler  shift  through  an  ionospheric  region 
which  varies  slowly  with  path  position  and  time  (Weeks,  1958;  Kelso,  1960,  1964; 
Ginzburg,  1964;  Bennett,  1967).  It  is  evident,  from  Eqs.  (46)  and  (49)  that  Eq.  (49) 
is  valid  only  when  the  Doppler  shift  is  small  compared  with  the  transmitter  fre¬ 
quency.  For  most  problems  of  propagation  through  the  earth's  ionosphere,  Eq.  (49) 
is  an  adequate  approximation  for  the  instantaneous  frequency.  However,  there  are 
laboratory  plasmas  and  some  planetary  atmospheres  (for  example.  Jupiter)  where 
(49)  may  not  be  a  good  approximation.  In  addition,  the  constraint  of  Eq.  (46)  may 
not  hold  in  many  space-time  varying  dielectrics.  In  the  next  section,  we  will  study 
the  exact  solutions  of  Eqs.  (11)  and  (14)  in  some  dielectric  materials. 
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Equation  (52)  is  not  simply  solved  for  arbitrary  variations  in  <  with  z  and  t.  How¬ 
ever,  for  specific  variations  in  *  solutions  are  readily  obtained. 

3.1.1  0  SEPARABLE 

When  «  is  separable  we  may  write  0  =  0^  (z)  /32(t).  In  this  case,  upon  com¬ 
bining  the  first  and  third  members  in  (48)  we  get 

w/32(t)  =  C1.  (53) 

Therefore,  the  observer  moving  with  the  velocity  V  =  1/0  finds  that  w/3^  (t)  is  an 
invariant.  Similarly  combining  the  first  and  second  members  of  (52)  we  get 

fZ  0,  (z')  dz!  -  f4 — •  =  C„  ,  (54) 

J  1  J  jS2(f)  2 

so  that  the  general  solution  of  (51)  when  <  is  separable  is 

^  )d«'  -/—]  •  » 

2 

where  the  arbitrary  function  4>  is  determined  by  specifying  boundary  conditions  on 
(o  along  any  curve  in  the  z  -  t  plane.  For  example,  if  0 2  did  not  depend  on  time, 
and  we  specified  w  =  a>0  at  z  =  0  for  all  t,  then  <t>  =  constant,  and  therefore 
tu  =  constant,  as  would  be  expected. 

3. 1.  2  TAYLOR  EXPANSION  OF  0 

The  situation  when  0  is  separable,  studied  above,  does  not  usually  occur  in 
practical  situations.  However,  there  are  often  problems  in  which  it  is  appropriate 
to  Taylor  expand  0  (z»  t)  in  either  z  or  t.  For  example,  we  can  consider  the  case 
in  which  we  desire  to  study  the  propagation  only  over  the  time  interval  tj  £  t  <  tg. 
In  that  interval  we  may  expand  0  in  Taylor  series  in  t  as 

0  (z,  t)  =  0Q  (z)  +  t  0  j  (z)  +  •  •  •  .  (56) 

If  we  assume  that  the  first  two  terms  of  the  Taylor  aeries  are  an  accurate  repre¬ 
sentation  of  0  in  the  interval  tj  S  t  <  tg  we  obtain  from  Eq.  (52) 

dz  dt  dm 

1  '  0(z,  t)  =  '  w01  (iz)  * 


(57) 


Prom the  first  and  third  members  o'.  Eq.  (5?)  we  obtain 


S(z)  r, 

*  C« 


St 

v'havs  S{z)  =  j  /5j(s' )  dz1  .  Therefore  <u  exp[S(z)]  is  an  invariant  as  one  moves 
with  the  group  velocity  \i!i  .  A  second  invariant  is  obtained  by  solving  the  first  and 
second  members  of  f 57).  ThisgU*'? 


t»'S<8)  -  /* dz'  0q(t.'  }  e“S{s' }  «  C?  . 

Upon  using  Eqs.  (58)  and  {59)  in  the  general  solution  Cj  =  $ (C2)  wa  get 

»(*.  t)  .  e-Sl8>4  -/V  ^!z')e-s«s'»]  . 


The  wavenumber  x(a,  t)  is  related  to  u  by  k  =  u|3  .  We  can  determine  the  arbi¬ 
trary  function  <?>{•••)  by  specifying  boundary  conditions  on  w .  For  example,  sup¬ 
pose  we  specify  that  w  =  ioQ,  for  all  t,  at  z  =  0.  (This  condition  is  appropriate  for 
the  case  of  a  plane  wave  of  frequency  uQ  transmitted  into  a  space -time  varying 
half-space. )  This  requires  that,  in  Eq.  (56),  the  function  4>  =  constant.  We  there¬ 
fore  obtain 


u  -  wQ  e 


-J  £,(z')dz' 


In  the  limiting  case  when  (3Q  and  /3j  are  independent  of  z  this  result  can  be  readily 
shown  to  be  identical, with  the  previous  result  of  Morganthaler  (1958)  that 


“  ■  “o  0(t)  '  “o  (  Ho  ;  0tt)  ■ 


where  tQ  is  the  time  at  which  the  signal  at  (z,  t)  was  at  z  =  0.  To  see  this,  we 
recall  that  in  spatially  homogeneous  dielectrics 


f  dt*  f  dt1 


mrff  rtratmw<ra©9*33Sl ! 
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Upon  performing  the  integral  in  (63a)  we  obtain 


"W«wv 


(63b) 


if  (63b)  is  substituted  into  (62)  it  is  evident  that  Eq.  (61)  and  (62)  are  identical  (in 
the  limit  of  spatially  homogeneous  dielectrics).  It  is  also  interesting  to  use  (61) 
in  Eq.  (38)  to  calculate  the  velocity  V^)  with  which  values  of  k  are  propagated  whan 
/3  is  given  by  Eq,  (56).  Using  (61)  in  (38)  we  find 


A) 


.d. 


•i 


[**i  ~($\  ] 


(64) 


3. 1.  3  TRAVELING- WAVE-MODULATED  DIELECTRICS 

One  other  case  of  interest  is  when  $  depends  on  position  and  time  as  z  -  vQt. 


This  situation  arises  when  the  dielectric  is  modulated  by  a  travelling  wave,  such  as 
an  acoustic  wave.  For  th 
two  members  of  Eq.  (52) 


an  acoustic  wave.  For  this  case,  upon  defining  u  =  z  -  vQt,  we  have  from  the  first 


dz  yt  dz  -  vodt  du _ 

1  _  vq/3  '  1  -  vQP  '  1  -  vQp 


(65) 


Combining  the  last  member  of  (65)  with  the  third  member  of  Eq.  (52)  gives 


du 


1  -  v  5(u) 


dw 


d  B 

V  <o  i — 
o  du 


(66) 


Equation  (66)  is  readily  solved  to  give 


( 1  -  (z  -  vQt)  ]  ui(z,  t)  =  Cj  . 


(67) 


Therefore  w(l  -  vQ0)  remains  invariant  as  one  travels  with  the  group  velocity 
V  =  1//3  .  It  is  also  simple  to  obtain  the  other  constant  of  motion  Cg  which  is 
determined  from 


dt_ 

dz 


«  0  <*  -  vQt)  . 


(68) 


Letting  ->  r  .  v +  eo  „  jit  .  _gu_ 

0  o  dz  1  dz  ‘  we  flnd  as  the  solution  to  (68) 


(z  - vot) 


!■ 


_,dS _ 

v0/3(S)  “  z  =  S  * 


The  general  solution  for  «  is  ther  efore 
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m 


$> 


<"'z.  t> 


Z  -  V  t 

[/( rk-. 


»)- 


vo/2iS) 


(70) 


z::t:z tz?  * ,s 

*•  ft  -  ITT- ‘JZT'tTT* "ave  distupb"eB 

performing  the  Integration  (a«8umI„g  jj  .i  |  ,  s”\  j  <™>  "'er 


<5 


w  (z,  t)  = 


/Tv 


Jin  a(z  -  vQt) 


tun 


fci 


1  *  7~*  1  + 


Si — f 

—  COS  a  (2  -  v  t) 
P  /  ° 


(71) 


if  we  again  specify  that  u  '  »  at  a  »  o  for  ,  u 

LT  8toMi- ,fter  “n“  1;:."^;.,:;%;^’ 


<j>(p)  =  w 


'  vp)(* 1  r,i)*  (-^jyTTrp2 


■<p2;>  „ 

p  1 


v-a) 

V  j 


(72) 


where 


and  p  is  the  argument  of  $('  •  • )  in  Eq.  (71).  It  appears  from  Eqs.  (71)  to  (73)  that 
w  (z,  t)  can  exhibit  some  rather  unusual  behavior  in  the  limit  when  vQ  —  Vp .  This 
problem  has  been  considered  by  He-esel  and  Oliner  (1961)  who  have  shown  that  the 
difficulty  arises  because  v  —  v  gives  rise  to  a  singularity  in  the  basic  differential 

"  r 

equation  of  Eq.  (la).  That  is,  if  (la)  is  written  in  terms  of  the  variable  s  =  z  -  vt, 

then  the  coefficient  of  the  dElda  term  becomes  vanishingly  small  as  vQ  —  Vp. 

Therefore  our  results  in  (71)  to  (73)  are  limited  to  the  regions  where  vQ  is  not  too 

close  to  v. 

P 

3.2  Time  of  Transit  in  Space-Time  Varying  Dielectrics 

Sitvte  the  phase  and  group  velocities  in  a  space -time  varying  medium  are  func¬ 
tions  opposition  and  time,  it  is  not  immediately  evident  how  long  it  would  take  for 
a  disturbance  to  travel  a  distance  L  along  a  ray.  To  study  this  problem  let  us  con¬ 
sider  the  motion  of  the  point  at  which  the  phase  =  0.  In  particular,  let  us  suppose 
<!>  =  0  for  t  *  tQ  at  some  point  ^  on  a  given  ray.  Then  the  time  tL  at  which  4>  =  0 
will  reach  .mother  point  along  the  ray  is  a  solution  of 


($L*  *L} 

J (k  •  dx  -  w 


dt)  =  0 


«5o-  V 


(74) 


To  study  the  solution  of  Eq.  (74)  let  U3  specialize  to  the  case  when  k  is  in  the  same 
direction  as  V«  .  If  we  assume  that  the  initial  point  x^  is  z  =  0  and  the  observation 
point  is  at  z  =  L  we  can  rewrite  Eq.  (74)  as 


(75) 


Since  the  line  integral  in  Eq.  (75)  is  independent  of  path  it  can  be  taken  along  any 
curve  joining  (0,  tQ)  to  (L,  t^)  in  the  z-t  plane.  For  instance,  we  shall  find  it  con¬ 
venient  to  write 
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^k(z,  o)  dz  -  f  w(L,  t)  dt  =  0  . 


4> 


(76) 


Equation  (76)  is  an  integral  equation  to  be  solved  for  the  transit  time  (t,  -  t  ),  To 

Jj  o 

understand  the  meaning  of  Eq.  (76)  let  us  first  suppose  that  /3(z,  t)  is  independent 
of  t,  and  that  w  (z  =  0)  is  specified  to  be  equal  to  wQ  for  all  t.  Then  since 
<o  (z  =  0)  =  w(z  =  L)  =  «0,  Eq.  (76)  become**  .he  obvious  result 


4-1 

-o -!■ 


dz 


V(z)  ' 


(77) 


-1/2 

where  V(z)  =  [  fiQ<  (z)]  .  Ir  (he  other  limit  when  {}  (z,  t)  is  independent  of  z.  we 

get  (assuming  k(t  =  t  )  =  k  for  nil  zl  that  tT  is  a  solution  of 

O  O  L 


,l 

«  f  V(t’)dt 


(78) 


which  is  the  reuvlt  obtained  previously  (Pante,  1971). 

As  an  example  of  the  application  of  Eq.  (75)  to  dielectrics  which  vary  with  both 
position  ar  i  time,  let  us  consider  the  case  when  /3  is  given  by  Eq.  (56),  along  with 
the  boundary  condition  that  w  =  wQ  at  z  =  0  for  all  t.  Then  from  Eq.  (61)  we  have 

-  J  /3j(z' )  dz' 

w  =  wQ  fe  0  ,  (79) 


k  -  /3(z,  t)  e 

V* 


"/  0i<z')dz' 


(80) 


Using  these  equations  in  (75)  gives  for  the  transit  time 


A..tL-t0.eS«'>/V,>e'S(*',dz'. 


(81) 


where  S(z)  =  /?j(g)dg.  From  Eq.  (81)  one  could  also  define  an  equivalent 

velocity  V(L)  -  L/At,  as 


group 
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V  =  - ; - - - -  .  (82) 

eS(L)  ^0(Z')e-S(z,)',z 
Jo 

In  the  limit  of  spatially  homogeneous,  tin:^- invariant  dielectrics  ( /?.  =0,  /?  inde- 

-  -1  *  ° 

pendent  of  z)  Eq.  (82)  reduces  to  the  usua.’  result  V  •-  (,Q)  . 


3.3  Discussion  of  Transmission  Through  a  Dielectric  Slab 

The  results  of  Sections  3. 1  and  3.  2  can  be  app.  to  considering  the  transmis¬ 
sion  of  plane  waves  through  a  lossless  dielectric  slab  of  thickness  L  in  which  the 
permittivity  varies  slowly  with  space  and  time.  Let  us  denote  the  solutions  in  the 
slab  by 


E  = 


eo(z,  t; 


V 


/( wdt  +  kdz) 


(83' 


where  w  and  k  are  the  appropriate  solutions  of  Eq.  (ST)  and  eQ(z,  t)  is  the  appropri¬ 
ate  solution  of  Eq.  (39).  For  example,  if  /3(z,  t)  =  f) j(z)  /32<t)  then 


eo(Zt  t;  V  to)  =  eo(0) 


-*(./ 

1/2 

fWl 

_73jlS> 

[A<»  J 

3/2 


(84) 


Now  suppose  the  slab  occupies  the  region  OS  z  S  L,  and  a  plane  wave 
Ej  =  exp  (iwQ(t  -  z/c)]  is  normally  incident  upon  the  slab  from  the  region  z  <  0. 
Then  the  field  transmitted  through  the  slab  at  time  tg,  into  the  region  z  >  L,  v/ill 
consist  of  a  number  of  components.  First,  there  will  be  a  wave  which  ci  osses 
(from  z  <  0)  the  z  =  0  boundary  at  time  t  j,  and  is  transmitted  directly  through  the 
dielectric,  arriving  at  z  -  L+  at  the  time  t  .  Next,  there  is  a  wave  which  crosses 
the  z  =  0  boundary  at  time  r.,.  and  arrives  at  z  =  L  at  the  time  t  ,  after  being 
internally  reflected  at  time  Tg  by  the  z  =  L  boundary  and  at  time  Tj  by  the  z  =  0 
boundary.  Next,  there  is  a  wave  which  crosses  (from  z  <  0)  the  z  =  0  boundary  at 
time  r,  and  arrives  at  z  =  L+  at  the  time  t  ,  after  being  reflected  twice  at  the  z  =  0 
boundary  (at  times  ^  and  r?)  and  twice  by  the  z  =  L  boundary  (at  times  v 2  and  r^), 
and  so  forth.  Let  us  define  R(0,  r  )  as  the  internal  reflection  coefficient  at  z  -  0+ 
boundary  at  rime  r,  K(L,  r  )  as  the  internal  reflection  coefficient  at  z  =  L’  at  time 
t  (for  the  case  when  the  medium  is  spatially  homogeneous  K(0,  r )  =  R(L,  r  ) 

=  [M7\lto  -  lj  [\Mr )/  «0  +  l]’1),  T(r)as  the  transmission  coefficient,  from 
z  £  O'  to  z  =  04  -t  time  r,  and  T(r)  as  the  transmission  coefficient  from  z  =  L'  to 
u  -  L4 ,  Then,  the  transmitted  field  at  z  =  L4  can  be  written  as 
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E(z  -  L+,  t  =  ts)  = 


T(ts) 


lb) 


T(Tj)  A(Tj)  e 


o  1 


+  T( Tg)  A(t3)  R(0,  Tj)  R(L,  r2)  e1  “ 


ot3 


(85) 


+  T(r5)  A(t5)  R(0,  Tj)  R(L.  t2)  R<0,  Tg)  R(L,  rj  e  "°T5  +  ...}, 

where  A(r)  =  eQ(z  =  L+,  t  =  tg;  zq  =  0,  tQ  =  t).  From  (85)  it  is  clear  that  the  nature 
of  the  transmitted  field  will  be  known  once  the  times  r^,  t2>  Tg  •  •  • ,  have  been 
determined.  Extending  the  discussion  of  Section  3.  2  we  see  that  these  are  solutions 
of  the  equations 

L 

J  k(z,  tp)  dz  = 
o 

o 

-  J k(z,  i-j)  dz  = 

L 

L 

J k(z,  Tg)  dz  - 
o 

Therefore  Tj,  r2>  and  so  forth,  could  be  obtained  by  first  solving  (86)  for  rr 
Then  using  the  solution  of  Eq,  (86)  in  Eq.  (87),  the  resulting  equation  could  be 
solved  for  Tg,  and  so  forth. 

Equation  (85>  can  also  be  applied  to  the  case  of  transmission  into  an  infinite 
half-space.  We  can  obtain  this  limit  by  setting  T(t)  =  1  and  R(L,  r  )  =  0.  We  there¬ 
fore  obtain  for  the  field  at  a  point  z  at  time  t 

s  s 

E<V  V  =  T(VA<Ve  °  1  '  (89) 


where  Tj  is  the  solution  of  the  equation 
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w<0,  t')  dt’  . 


(90) 


When  /9(z,  t)  is  given  by  Eq.  (56),  wt  'an  use  the  results  of  Section  3.  2  to  express 
as 


(91) 


Once  Tj  is  known,  T(tj)  can  be  obtained  by  applying  the  usual  WKB  methods  to  the 
profile  ( (z)  =  /iQ'1  [  0Q(z)  +  i-j  p^z)]2. 

As  an  example  of  the  application  of  the  more  general  result  of  Eq.  (85)  let  us 
consider  the  case  when,  the  dielectric  slab  is  spatially  homogeneous.  Then  (85) 
becomes 

A  r  i«  r. 

E(z  =  L+.  t  =  tg)  =  T(ts)  |t(tj)  A<y  e  °  1 


+  T(t3)  A(t3)  R(Tj)R(T2)e 


o  3 


(92) 


+  T<  y  A(  r5)  R(  r2 )  R(  r2)  R(  y  R(  r4)  e  °  5  +  “  •  •  }  , 


where 


A(r) 


3/4 


and  r  is  the  solution  of 
m 


v(t' )  dt' 


mL 

m  =  1,  2,  3,  •  •  • 


(93) 


(94) 


To  examine  the  various  frequency  components  present  in  the  transmitted  wave  of 
Eq.  (92)  we  can  Taylor  expand  the  functions  rm(t).  That  is 
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(95) 


(96) 


Using  this  result  we  may  rewrite  (92)  as: 


E(z  =  L+,  t) 


(t  -  ts> 


(97) 


where 


Bl  - 


<(r21-l) 


3/4 


i  w  r . ,  , 
o  21-1 


2(1-1) 

n 

j=0 


R(r.) 


and 


'  v<y 

“0 

1 

< 

to 

l 

“0 

L  €(v  J 

(98) 


From  Eq.  (97)  we  see  that  the  transmitted  signal  consists  of  components  at  the 
instantaneous  frequencies  ftj,  fig,  fig,  Sl^"’  •  This  interpretation  assumes 
that  B, j  is  a  slowly  varying  function  of  time  in  comparison  with  exp  (ift^t),  which 
is  the  case  in  the  WKB  approximation.  The  importance  of  the  frequency  cc  oonents 
at  Sly  ^3»  ^4*  and  so  forth,  in  comparison  with  that  at  fi ^ ,  will  depend  on  the 
amplitude  of  the  reflection  coefficient  B.  For  R  «  i  only  ft,  will  be  significant, 
but  for  R  near  unity  this  conclusion  is  clearly  not  true. 

Therefore,  we  see  that,  because  of  the  spatial  boundaries,  the  transmitted 
signal  has  components  at  ft  ,  ftg*  ^3  "  ’ '  *  and  no*  Just  a*  as  found  by 
Morganthaler  (1958),  who  did  not  account  for  boundary  effects.  For  relative  per¬ 
mittivities  near  unity  the  components  at  Sly  ftg,  and  so  forth,  will  be  negligible 
compared  with  that  at  ftj.  For  large  relative  permittivities,  however,  the  higher 
order  frequency  components  will  be  significant. 
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3.4  Further  Discussion  -  Day  Tracing  Methods 

In  the  general  case  when  the  direction  of  propagation  is  not  along  V< ,  it  is 
not  a  simple  matter  to  solve  Eq.  (11)  or  (14)  for  space-time  varying  dielectrics. 

In  many  cases  it  is  acceptable  to  use  the  approximation  of  Eq,  (49),  but  for  others 
this  may  not  be  possible.  In  such  cases  it  appears  most  appropriate  to  approxi¬ 
mate  the  temporal  behavior  of  the  dielectric  in  a  stepwise  fashion.  For  example, 
in  the  time  interval  0  £  t  £  t^  the  function  (s  (x,  t)  can  be  approximated  by: 

P  =  p1  (x)  for  0  £  t  <  tj  /3  =  (x)  for  tj  <  t  <  t2  •  •  •  ,  &  =  £N  (x)  for  tN_  2  <  t  <  tN> 

where  we  assume  that  |  { >9j  -  (3^  ^)!  (2^  |  <<  1.  Now  suppose  we  have  a  signal  of 
frequency  Wj  which  enters  the  dielectric  along  the  direction  kj  at  time  t  =  0.  Then 
in  the  interval  0  £  t  <  tj  we  can  ray  trace  (using  Snell's  Law)  to  determine  the  path 
along  which  the  signal  will  travel.  The  distance  Sj  along  the  ray,  which  the  phase 
<f>  =  </>o  will  traverse  in  time  interval  tj,  is  obtained  by  solving 


t 


1 


0j(x)  dS  , 


(99) 


where  dS  =  k  •  dx,  and  the  integral  is  taken  along  the  ray  path.  At  t  =  t,,  P  is 
suddenly  changed  to  /3g(x).  The  new  frequency  associated  with  the  point  we  are 
following  along  the  ray  is  (see  Eq.  33) 


/Ux.) 

<J  -  (i)  -  -  ■  +  ■■  *■  .■ 

2  ’  1  P2(xx)  ’ 


(100) 


where  Xj  is  the  coordinate  of  the  point  S  =  Sj  on  the  ray.  Also,  if  {c(tj  -  8)  is  the 
normal  to  the  phase  surface  at  tj  -  8  ,  then  at  tj  +  8  the  normal  is  still  in  the  same 
direction,  by  virtue  of  Eq,  (31).  That  is,  if  the  angle,  between  k(tj  -  8  )  and  the 
z  axis  is  9 j,  then  the  angle  between  k{tj  +  8  )  and  the  z  axis  is  also  8 j.  The  new 
ray  path  (valid  during  t.  <  <  <  t„)  *s  deteri.v.naa  by  ray  tracing  starting  at  x  =  x., 

*  W  A  -1 

along  the  initial  direction  k(tj  t-  8  )  =  k(tj  -  8  ).  The  distance  S2  -  Sj  which  the 
wa»efront  wilt  travel  along  the  new  ray  path  in  the  interval  tj  <  t  <  tg  is  obtained 
by  solving 


t 


2 


(101) 


We  can  continue  the  above  process  until  we  have  traced  the  progress  oi  the  ray 
point  we  have  chosen  for  all  values  of  time  in  the  interval  0  S  t  £  . 
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4.  APPLICATION  TO  DISPERSIVE  MEDIA 


4. 1  Discussion  of  the  Invariants 


The  propagation  in  dispersive  media  is  far  more  complex  than  in  dielectrics. 
As  the  simplest  example,  consider  the  propagation  in  a  lossless,  isotropic  plasma 
in  which  the  propagation  is  in  the  same  direction  as  Vu  ,  where  to  is  the  elec- 

p  p 

tron  plasma  frequency.  In  this  case  the  dispersion  relation  is  c2k2  =  w2  -  w  2(x,  t), 
so  that  Eq.  (14)  becomes  (assuming  to  is  not  too  close  to  wp)  P 


dto .  1  duip2 

7z~  KT  “3t  • 


(102) 


From  Eqs.  (17)  and  (18)  we  have  that  the  solution  of  (102)  will  have  the  form 
Cj  =  ^(Cg),  where  Cj  =  f(w,  z,  t)  and  Cg  =  g(  to,  z,  t)  are  the  particular  integrals  of 


dt  _  dz _ _  2w  dw 


(103) 


Unfortunately,  Eq.  (103)  is  not  readily  solved  for  the  constants  of  motion,  except 
in  the  limits  when  aip  depends  on  z  or  t,  but  not  on  both.  For  example  (as  expec¬ 
ted  from  the  discussion  of  Section  2)  we  have  that  to  =  Cj  when  to  is  independent 
of  time.  This  means  that  to  la  an  invariant  as  one  moves  with  the  group  velocity 

V  =  c(l-wp/w  )  .  When  to  is  independent  of  position  we  find  that  the  particular 

2  2 

integral  Cj  =  to  -  «p(t),  which  is  just  a  statement  of  the  fact  that  values  of  k  are 
propagated  with  the  group  velocity  on  spatially  homogeneous  media. 

We  have  not  been  able  to  obtain  the  exact  invariants  of  the  motion  when  w  is 
an  arbitrary  function  of  z  and  t.  In  most  practical  plasma  problems,  however^  we 
are  fortunate  that  the  inequality  of  Eq.  (46)  is  satisfied  so  that  to  can  be  approxi¬ 
mated  by  the  result  of  Eq.  (49).  One  can  easily  imagine  problems  where  Eq.  (46) 
is  not  satisfied,  however,  and  in  this  base  ray  tracing  methods  seem  most  appropri¬ 
ate.  In  the  next  section  we  will  discuss  how  ray  tracing  can  be  applied  to  space- 
time  varying  plasmas. 


4.2  Ray  Tracing  in  Space-Time  Varying  Dispersive  Media 

The  simplest  method  of  applying  ray  tracing  techniques  to  a  lossless,  isotropic, 
space-time  varying  plasma  is  to  approximate  the  temporal  variations  in  a  stepwise 
fashion  as  was  done  in  Section  3.  4.  That  is,  we  approximate  the  plasma  frequency 
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a)p(x,  t)  by:  Wp  =  "pl  W  for  0  <  t  <  tj,  wp  =  u>p2(x)  for  ^  <  t  <  tg,  tup  =  tup3{x)  for 
tg  <  t  <  tg,  and  bo  forth.  We  assume  that 


V  ~  v 

“pi 


«  1 , 


and  so  forth.  To  illustrate  the  method  we  will  further  assume  that  the  plasma 
occupies  the  half- space  z  >  0,  and  is  spatially  stratified  in  the  z  direction  only. 
At  t  =  0,  a  point  p  on  the  envelope  of  a  signal  of  frequency  Wj  enters  the  medium 
at  an  angle  relative  to  the  z  axis.  Then  by  Snell's  Law  we  have  that  the  ray 
path  during  the  time  interval  (0,  tj)  is  determined  from 


u  2  0) 
i  - pl ; 

X  /  & 

sin  8  - 
0 

", .  v'f  ’ 

1  2 

1  2 

L  wi  j 

L  wi  J 

sin  Q(z) 


The  component 


* 


at  frequency 


tdj  moves  along  this  path  at  the  group  velocity 


=  c 


|l  -  ui 


2,  2 

pi  '  W1 


11/2 


so  that  during  the  time  interval  (0,  tj)  the  distance  Sj  travelled  along  the  ray  by  the 
point  p  on  the  envelope  of  the  signal  is  the  solution  of 


Sj 

r_ds_ 

J  V,  (z)  * 


(104) 


where  the  integration  is  along  the  ray  path  (in  isotropic  media  the  ray  and  group 
paths  are  the  same).  We  could  also  have  chosen  to  follow  the  progress  of  a  par¬ 
ticular  value  of  phase  <fi  -  <p  j.  In  the  time  interval  tj,  the  distance  Sj '  travelled 
along  the  ray  by  this  value  of  phase  is  then  the  solution  of 


/ 


dS 

v. 


Pi 


(105) 


#The  above  discussion  would  appiy  to  the  motion  of  a  wave  packet  with  fre 
quency  spectrum  initially  centered  sharply  about  ui  = 
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where  vpJ  =  c/yi  -  Upl/  . 


Now  suppose  that  at  t  =  tj  -  8  ( 8  -*  0)  the  point  p  on  the  envelope  has  reached 
the  point  x  =  xJt  and  that  the  ray  path  at  that  point  makes  an  angle  8^  with  the 
z-axis.  At  t  =  tj  the  plasma  frequency  is  suddenly  changed  from  wpl(z)  to  wp2(z). 
By  virtue  of  Eq.  (31),  8  cannot  change  instantaneously  so  that  fi(tj  -8)=0(tj+8) 
=  By  Therefore,  the  new  ray  path  in  the  plasma  during  the  interval  t^  <  t  <  tg  is 
determined  from 


1  - 


°2  J 


1/2 


sin  =  j  1  - 


«p2<z> 


■^2  1/2 


2  J 


sin  C  (z)  , 


(106) 


where  the  new  frequency  associated  with  the  point  p  we  are  following  is 

1 1/2 


W2  = 


2  2,  ,  2,  v 

“l  -  V(Z1)+  "p2(zl) 


(107) 


This  component  now  moves  along  the  ray  trajectory  with  the  group  velocity 

,1/2 


v„ 

=  V„  =  c 

1  - 

—2 

2 

l  w2  /  J 

(108) 


and  during  the  interval  tj  <  t  <  t2  traverses  the  distance  S2  -  Sj  along  the  ray, 
where  S2  -  Sj  is  the  solution  of 


(109) 


and  the  integral  is  along  the  ray. 

The  above  procedure  can  be  repeated  continually  until  the  position  and  fre¬ 
quency  of  the  point  p  of  interest  have  been  determined  everywhere  in  the-time 
interval  of  interest.  Of  course,  if  in  any  region  of  space  (or  time)  we  reach  the 
situation  where  w  is  close  to  cup,  the  WKB  method  is  no  longer  valid,  and  one 
must  do  a  more  careful  analysis  (see,  for  example,  Kelso,  1964;  Ginzburg,  1964). 
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5.  COMMENT  ON  ABSORPTION 


To  keep  our  discussion  relatively  simple  we  have  avoided  considering  the 
effect  of  absorption.  When  absorption  is  present  we  find  that  [since  (x,  t)  must 
be  unique  even  when  losses  are  present]  Eq.  (9)  still  holds,  except  that  w  and  k 
may  both  be  complex  numbers  (and  in  addition  the  real  and  imaginary  parts  of  the 
k  vector  will  generally  lie  along  different  directions).  However,  our  interpretation 
of  various  quantities,  such  as  the  group  velocity  and  the  ray  paths  is  no  longer 
valid.  For  example,  when  absorption  is  present  |  V  |  =  j  |  may  be  greater 
than  the  speed  of  light,  and  consequently  no  longer  retains  any  physical  meaning 
even  in  the  limit  of  spatially-homogeneous,  time  invariant  media. 

In  order  to  understand  why  w  and  k  may  both  be  complex  in  lossy  space-time 
varying  media,  we  first  realize  that  since  Eq.  (9)  is  valid  when  Eq.  (31)  must  also 
still  hold,  so  that  if  the  properties  of  a  lossy  medium  are  abruptly  changed  at 

t  =  t.  we  still  have  that  k(t,  *-  8  )  =  k(t,  -  S  ).  In  a  lossy  dielectric  k  and  m  are 

*  2  2~  *  *  * 

related  through  k  -  m  nQ*  +  iw  /iQ<r  ,  where  a  is  the  conductivity  of  the  medium. 
We  now  suppose  that  «  =  «j,  <r  -  o-j,  and  w  =  Wj  =  real  for  t  <  tjj  and  «  =  «2, 
o-  -  o-g,  for  t  >  tj.  Then  since  k  =  j  k  j  is  continuous  at  t  =  tj  we  can  determine  the 
new  frequency  v<-  during  t  >  t.  from 

It  1 

2  2 

“1  M0*l  +iall^ocrl  =  w2  ^o«2  +  iw2^oa2  •  (110) 


If  u>2  were  real  we  would  require  that 


and 


(Ilia) 


(111b) 


Equations  (111a)  and  (111b)  cannot  be  simultaneously  satisfied,  except  for  the 
special  case  when  ( « ^ /  -  (o-j /  c r^).  We  therefore  conclude  that  for  Eq.  (110) 

to  be  satisfied  in  general,  mus*  be  complex.  In  particular  denoting 
gl>2  =  w2  *  we  ^ave 


tl 

w2  ®2  ’ 


(112a) 
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wl°l  =  2u>2 


«2  + 


2  2  • 


(112b) 


Equations  (112a)  and  (112b)  can  be  solved  for  Wg  and  w'g.  We  will  not  probe  too 
deeply  Into  the  effects  of  absorption.  This  will  be  deferred  to  a  later  paper.  We 
do  comment  that  the  analysis  becomes  relatively  straightforward  in  the  limit  when 
the  absorption  is  small  (cr/ut  «  1)  since  then  the  relevant  equations  may  be 
studied  using  perturbation  methods.  For  example,  the  ray  paths  can  be  determined 
by  neglecting  absorption,  and  so  forth. 


6.  SUMMARY 

We  have  studied  the  properties  of  the  WKB  solutions  in  lossless,  isotropic 
space-time  varying  media.  It  was  found  that,  in  principle,  one  can  always  obtain 
constants  of  the  motion  which  lead  to  a  complete  determination  of  the  frequency 
and  wavenumber,  once  appropriate  boundary  conditions  have  been  specified.  Once 
tu  and  k  are  known  one  can  readily  study  the  transmission  through  space-time 
varying  media,  as  was  illustrated  in  Section  3.  3.  In  general,  however,  the  con¬ 
stants  of  motion  are  not  always  easily  obtained.  In  such  problems  we  have  shown 
in  Sections  3.4  and  4-  2  that  one  can  obtain  solutions  by  modelling  the  medium  by  a 
series  of  temporal  s*eps.  That  is,  the  index  of  refraction  n(x,  t)  is  approximated 
by  n^x)  in  0  i  t  i  tj,  n2(x)  in  tj  <  t  <  t9,  and  so  forth,  and  ray  tracing  techniques 
are  applied  during  each  time  interval. 
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